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We discuss some properties of the ferromagnet - superconductor proximity sys- 
tem. In particular, the emphasis is put on the physics of the Fulde-Ferrell-Larkin- 
Ovchinnikov {FFLO) like state. In addition to Andreev reflections it features a 
number of unusual thermodynamic and transport properties, like: oscillatory be- 
havior of the pairing amplitude, density of states and superconducting transition 
temperature as a function of the ferromagnet thickness. Surprisingly, under certain 
conditions spontaneous spin polarized current is generated in the ground state of 
such a system. We provide some informations regarding experimental observations 
of this exotic state. 

I. INTRODUCTION 

When a normal non-magnetic metal is connected to a superconductor it acquires super- 
conducting properties, like non-zero pairing amplitude. This effect, known as the proximity 
effect Q|, has extensively been studied for almost half a century. It is rather well under- 
stood by now in terms of Andreev reflections [2j], according to which an impinging electron 
(with energy less than superconducting gap) on the normal metal (NM) / superconductor 
(SC) interface is reflected back as a hole and the Cooper pair is created in superconductor. 
From the point of view of Andreev reflections the proximity effect can be regarded as a 
non-zero density of the Andreev correlated electron - hole pairs on the normal metal side of 
the interface. 



When a normal metal is replaced by a ferromagnet (FM), another energy scale enters 
problem, namely the exchange splitting which is related to the spin polarization of the 
electrons. Such FM/SC hybrid structures are important from the scientific point of view, 
as they allow the study of the interplay between ferromagnetism and superconductivity [3] 
as well as of device applications in such areas of technology as magnetoelectronics jj| or 
quantum computing [5|. 

It is widely accepted that ferromagnetism and superconductivity are two antagonistic 
phenomena, so one could expect that the proximity effect in FM/SC system should be 
suppressed. Indeed, the one can argue that in ferromagnet there are different numbers of 
spin-up (majority) and spin-down (minority) nj_ conduction channels, and due to the fact 
that incident and reflected particles occupy different spin bands, only a fraction njn^ of 
majority particles can be Andreev reflected Q]. 

On the other hand if an exchange field acts on the Cooper pairs, one would expect 
that either it is too weak to break the pair, or it suppresses completely superconductivity. 
However when a Cooper pair is subjected to the exchange field, it acquires a finite momentum 
and for certain values of the exchange splitting a new superconducting state is realized, 
known as Fulde - Ferrell - Larkin - Ovchinnikov (FFLO) state [?], Isj]. Interestingly such 
state features a spatially dependent order parameter corresponding to the non-zero center 
of mass motion of the Cooper pairs. This state features in non-zero spin polarization, 
almost normal tunneling characteristics and almost normal Sommerfeld specific heat ratio, 
anisotropic electrodynamic properties. Unfortunately the bulk state is very sensitive to the 
impurities and shape of the Fermi surface. Another novel feature of this state is a current 
flowing in the ground state. The unpaired electrons tend to one portion of the 

Fermi surface so a quasiparticle current is produced. In order to satisfy the Bloch theorem: 
no current in the ground state, a supercurrent, generated by the nonzero value of the pairing 
momentum, flows in opposite direction, and the total current is zero. 

Similar oscillations of the pairing amplitude have been predicted j^-^^ in ferromag- 
net /superconductor proximity systems. It turns out that these oscillations are responsible 
for the oscillatorybehavior of the SC critical temperature T c , first experimentally observed 
by Wong et al. [3], and the density of states ^| as the thickness of the FM slab is varied. 
In fact, the oscillations of the T r in FM/SC multilayers can be also explained in terms of the 
effective 7r-junction behavior It was shown that at specific FM thickness the Josephson 



coupling between two SC layers can lead to a junction with an intrinsic phase (of the order 
parameter) difference 5ip = tc, which exhibits a higher T c than the ordi nary one (dip = 0). 
The 7r-junction effect has been originally proposed by Bulaevskii et al. [15] to arise in the 
tunnel barriers containing magnetic impurities. It was also suggested that the 7r-junction 
can be realized in high-T c superconducting weak links Q], where the SC order parameter 
changes its sign under 7r/2 rotation. This has tremendous consequences as it leads to many 



important effects [17|,|18j], like: the zero energy Andreev states, zero-bias conductance peaks, 
large Josephson current, time reversal symmetry breaking, paramagnetic Meissner effect and 
spontaneously generated currents. 

From the point of view of the present paper the important issue is the formation of the 
Andreev bound states in FM/SC proximity system. The Andreev states arise due to the 
fact that the quasiparticles of the ferromagnet participating in the Andreev reflections move 
along closed orbits. Such states have been first studied by de Gennes and Saint- James 
in the insulator /normal metal/superconductor (I/NM/SC) trilayer. The energies of these 
states are always smaller than the SC gap A and symmetrically positioned around the Fermi 
level. They strongly depend on the geometry of the system as well as on the properties of 
the interfaces. In high-T c (<i-wave) superconductors, these states can be shifted to zero 
energy, due to the specific form of the symmetry of the order parameter Q], thus indicating 
7r-junction behavior in the system. Naturally, such Andreev states can also arise in the 
I/FM/SC heterostructures. Moreover, it is possible to shift the energies of these state s by 
changing the exchange splitting, as was first demonstrated by Kuplevakhskii & Fal'ko [21]. 
In turn, by properly adjusting the exchange splitting the position of the Andreev bound 
states can be moved to the Fermi energy. The system under such circumstances behaves 
like that being in the 7r-junction phase as the spontaneous current is generated [22]. 

Some of our results have already been published 22j- 2J|. Here we wish to present a 
more detailed study of the FM/SC proximity system in terms of FFLO physics. In some 
situations the ground state of FM/SC structures has properties of both the FFLO and the 
7r-junction, leading to various interesting and unexpected phenomena. 

The paper is organized as follows: In Sec. |H] the simple model which allows for self- 
consistent description of the FM/SC heterostructure is introduced. In Sec. IHII the nature 
of the Andreev bound states in the ferromagnet is discussed. The spontaneously generated 
current and corresponding magnetic field in the ground state are studied in the Sec. In 



Sec. El show some transport properties of the system, in the Sec. EUwe compare our system 
to usual FFLO state, and finally, we conclude in Sec. IVI1 



II. MODEL AND THEORY 

To study the properties of FM/SC system we have adopted the 2D Hubbard model 
featuring the exchange splitting in the ferromagnet and an electron - electron attraction in 
superconductor. The Hamiltonian is: 
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H = Yl Uj + \^E ex a - fi) 5ij cf a c ja + -J2 u i n i<? n i-° (!) 
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where in the presence of a vector potential A(r), the hopping integral is given by t. 
—te %e ^ ¥ i A ^ dr f or nearest neighbor lattice sites, whose positions are r*j and r}, and zero 
otherwise. The exchange splitting E ex is only non-zero on the FM side, unlike as Ui (electron 
- electron attraction) being non-zero only in SC. fi is the chemical potential, cf a , (cj CT ) are 
the usual electron creation (annihilation) operators and hi a = cf a Ci a . 

In the following we shall work within Spin - Polarized - Hartree - Fock - Gorkov (SPHFG) 
approximation [22j assuming periodicity in the direction parallel to the interface while work- 
ing in a real space in the direction perpendicular. Labeling the layers by integer n and m 
at each k y point of the Brillouin zone we shall solve the following SPHFG equation: 

C'^^JwA) = S nm S a p (2) 

where the only non-zero elements are: H^ m and = (lj — \oE ex ± // ± tcos(k y =)= 

eA(n)))S nm ±t<5„ in+ i for the upper and lower sign respectively, = and = FF^ m 
with a replaced by —a and H^ m = E^ m = —H^ = -H^ m = A n 5 nm and G"f, is corre- 
sponding Green's function (GF). 

As usual, the self-consistency is assured by the relations determining the FM (m n ) and 
SC (A n ) order parameters, current (J y i(i)(n)) and the vector potential (A y (n)) respectively: 

r, 2JV-1 

m n = n nT - n nl = - £ £ Re k y ) - Gf>„, k y ))e^ +1 ^ 2N } (3) 
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A n = U n J2(c nl (k y )c n ,(k y )) = yEE Re{GZ(^,k y )e^ +1 ^ 2N } (4) 



ky u=0 



Apt 2N ' 1 

J m) {n) = ^ £ sin(k y - eA y {n)) £ Re {G^ 33) K, k y )e^+^ N } (5) 

^ k y u=0 

A y (n + 1) - 2A y (n) + A y (n - 1) = -4vrJ^(n) (6) 
The details of the calculations can be found in 

III. ANDREEV BOUND STATES 

Before we discuss results of fully self-consistent calculations we would like to turn the 
attention to origin of Andreev bound states and take a look at physics of them from the 
point of view of semiclassical approach. 

From quasiclassical considerations, each bound state corresponds to quasiparticle moving 
along a family of closed trajectories [25j|. The energy of such bound state is determined by 
the Bohr- Sommerf eld quantization rules, according to which the total phase accumulated 
during one cycle has to be equal to multiples of 2tt. Interestingly, the bound states also 
emerge in the normal metal/superconductor (NM / SC) structures due to the Andreev 
reflections 3], according to which an incident electron is reflected back as a hole at the 
interface, and a Cooper pair is created in SC. Such states are built up from a combination 
of electron and hole wave functions. The example of the closed quasiparticle trajectory, 
producing the bound state, in an insulator/ (normal metal) /superconductor I /NM/ SC, is 
shown in the Fig. ^ It consists of an electron e segment, which includes a ordinary reflection 
at the I /NM interface, and hole h one, retracing backwards the electron trajectory. The 
total accumulated phase in this case consists of contribution from Andreev reflections at 
point A: —a\ + ipi and B: —a 2 + y?2 as well as contribution from the propagation through 
the normal metal (3(E). ai( 2 ) — arccos(E/\ A |) is the Andreev reflection phase shift, while 
<fi(2) is the phase of the SC order parameter at point A (B). /3(E) = 2L(k e — kh) + (3q is 
the electron-hole dephasing factor and describes the phase acquired during the propagation 
through the normal region, where the first term corresponds to the ballistic motion and the 
second one to the reflection at the I /NM surface. L is the thickness of NM, and k e (kh) is 
the electron (hole) wave vector. Thus the Bohr-Sommerfeld quantization condition is: 



- (at + a 2 ) ± (</?i - (p 2 ) + 13(E) = 2nvr 



(7) 




FIG. 1: The example of the quasiparticle path corresponding to the Andreev reflections, giving a 
bound state. The quasiparticle is trapped in the normal region because of normal reflection at the 
I /NM surface and the Andreev reflection at the NM / SC interface. The total phase accumulated 
during one cycle is equal: — (a\ + 02) ± {ip\ — ^2) + (3(E). 

where the ±(y?i — ^2) stands for the trajectories in the ±k y (parallel to the interface) 
direction. 

If there is no phase difference between points A and B in the Fig. ^(for example NM / SC 
interface), the bound states always appear in pairs symmetrically positioned around the 
Fermi level because of the time reversal symmetry in the problem. Moreover, due to the fact 
that there is no difference between electrons and holes at the Fermi level {(3{E = 0) = 0), 
there is no E = solution. In other words, the bound states always emerge at finite energies. 

The situation is quite different if there is a phase difference (ipi — between points A 
and B (see Fig. Q). The example can be the interfaces with c?-wave superconductors oriented 
in the (110) direction, where ((pi — ^2) = 7T- in this case, due to the additional phase shift ir 
the bound states can emerge even at zero energy. Such zero-ener gy A ndreev bound states, in 
the case of high-T c superconductors, have been predicted by Hu [20] and are known as zero- 
energy mid-gap states. The presence of the Andreev bound states at zero energy features in 
many important effects, like zero-bias conductance peaks, 7r-junction behavior, anomalous 
temperature dependence of the critical Josephson current, paramagnetic Meissner effect, 



time reversal symmetry breaking and spontaneous interface currents \v\ . Q| . 

Although the zero-energy states (ZES) are likely to appear when the phase of the order 
parameter at the interface is not constant, the resulting density of states at the Fermi energy 
is energetically unfavorable and any mechanism able to split these states will lower the energy 
of the system 3|- O n °f these is the self-induced Doppler shift 17, 27. | S = evpA, where 



A is a vector potential. The situation is schematically depicted in the Fig. EJ At low 
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FIG. 2: Generating of the spontaneous currents. 



temperature (T* « (£ /A)T c , where A is the penetration depth of the magnetic field) the 
splitting of the zero energy states produces a surface current. This current generates a 
magnetic field (screened by a supercurrent), which further splits ZES due to the Doppler 
shift effect. The effect saturates when the magnetic energy is equal to the energy of the 
Doppler shifted ZES. 
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lv the Andreev bound states also arises in I / FM / SC heterostructures 
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_ . _ ,., .32|- More importantly, as it was first predicted by Kuplevakhskii & Fal'ko 

^\ , it is possible to shift these states to zero energy by tuning the exchange splitting. So the 
crossing of the zero energy solution can be obtained either by changing the phase difference 
((/?! — (p 2 ) or by varying FM coherence length (exchange field). 

The properties of such bound states have been also studied fully quantum-mechanically 
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3l| and similar their behavior have 



within lattice models of the FM /SC systems 
)een obtained. Interestingly, it turns out, that as in the case of the high-T c structures 
such zero energy Andreev states support spontaneous currents flowing in the ground 



state of the FM/SC system 
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24]. The mechanism of generating of such currents 



is the same, as earlier discussed, namely the self-induced Doppler shift. So in fact, when 



the current flows, such one of the states will be twice shifted: once due to the exchange 
(Zeeman) splitting, and the second time due to the Doppler shift. 

For energies less than superconducting gap, the only Andreev bound states will contribute 
to the density of states p(E). However, as it was mentioned, for fixed thickness and exchange 
splitting, there will be Andreev bound states at different energies, for different angles of 
particle incidence (72 in the Fig. |TJ). Thus to get the density of states, one has to sum the 
energies of these states over all values of 72: 



tt/2 



(8) 



p(E)= ]T 5(E - E bound ) 

72=-V2 

and talk, in fact, about Andreev bands rather that single states. However, all that was 
said on properties of the bound states, remains true for Andreev bands too. In partic- 
ular the splitting of the whole band due to the spontaneous current is illustrated in the 
Fig. The additional structure comes from the other (higher order) Andreev reflections. 




FIG. 3: Doppler splitting of the zero-energy state. From Ref. 23]. 



Superconducting energy gap A = 0.376 in this figure. 

There is also a strong correlation between Andreev bound states (bands) and the pairing 



amplitude (22 



23, 



Each time the pairing amplitude at the I/FM interface changes its 
sign, the Andreev bound state (band) crosses the Fermi energy. Moreover in this case the 
spontaneous current is generated. 

From the experimental point of view the density of states, in particular its temperature 
dependence, can be a good measure of the current carrying ground state. At certain thick- 
nesses of FM for which the current flows there is a huge drop in the ptot( £ F) a t characteristic 



temperature T* ps (£s/A)T c , where £s and A are coherence length and penetration depth 
respectively. T* simply indicates the temperature at which magnetic instability, which leads 
to the generation of the current, takes the place. Such behavior is depicted in the Fig. |3]and 
should be observable experimentally. If there is no current the DOS is due to the Andreev 
band and is almost constant (we are well below T c ), and as soon as the current starts to 
flow the Andreev band splits so we observe a drop in ptati^F)- The important point is that 
T* and T c are different temperatures. 
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FIG. 4: The temperature dependence of the surface (FM/vacuum) density of states at the Fermi 
energy for various thicknesses of the FM slab in the figure. The solid (dashed) line corresponds to 
the solution without (with) the current. From Ref. [241 ] 




IV. SPONTANEOUS CURRENT 



The most remarkable feature of our calculations is that the solution of the SPHFG 
equations frequently converges to a solution with the finite current j y (n) even though the 
external vector potential is zero. The typical example of such a current, flowing parallel to 



the FM/SC interface, (j 



tot. 



n 



j y ]{n) + jyi(n)) is shown in the Fig. El for a few values 



of the exchange splitting. Behavior of the current, as a function of the layer index, is very 
similar to the density of states at the Fermi level. The oscillating nature of the current comes 



from the Friedel like oscillations of the DOS 
the DOS at the Fermi leve 
the current is very smooth 



23j . This is because current is proportional to 
. Within semiclassical calculations, which neglect these effects 
3J. 
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FIG. 5: The total (spontaneous) current jy 0t (n) = j y -\(n) + i v \ (n) flowing parallel to the FM/SC 
interface for a number of exchange splittings. From Ref. [2 j |. 

Another important issue is the distribution of the current through the whole trilayer 
structure. We find that it flows mostly in the positive y direction on ferromagnetic side and 
in the negative direction in the superconductor. Notably the total current, integrated over 
the whole sample, is equal to zero within numerical accuracy. This is as it should be for the 
true ground state and found to be in the FFLO state, where the current associated with 
the unpaired electrons is balanced by the supercurrent flowing in the opposite direction. 
Similarly here (see Fig. EJ). 




FIG. 6: Schematic view of the current distribution. 



Obviously, the spontaneous current distribution (see Fig. |BJ) generates the magnetic field 
through the sample. The total magnetic flux weakly depends on the thickness of the sample 



and the exchange splitting. Its magnitude is found to be a fraction of the flux quantum 
$o = h/2e and is smaller than upper critical field of the bulk superconductor. This is rather 
a large field and could be observable in temperature dependent measurements (see Fig. |7j). 

0.6 
0.4 



0.2 
0.0 

0.0 0.1 0.2 0.3 0.4 

T/T c 

FIG. 7: The temperature dependence of the total magnetic flux for thickness of the FM slab 
L/^s = 2.6 (solid), 6 (dashed) and 15 (dotted curve). 




V. TRANSPORT PROPERTIES 



Some information on spontaneous currents can be also obtained from conductance calcu- 
lations. To do so we attached a normal metal electrode to our FM/SC system and calculate 
current through NM / FM / SC system using nonequilibrium Keldysh Green's function tech- 
. To calculate t his current in terms of various physical processes we went along 



mque 



the way outlined in Ref. 



f. Q and got corresponding spin polarized formula for the current 
as a sum of four different contributions I — I\ + 12 + 13 + I a, where: 

h = ^ 2 t\ F e - J2jdu\l + G%M\ 2 PNNMPF F M[f{" - eV) - /( W )] (9) 

h = 4tt 2 4 f |E / dujlG^i^p^^pf^M^ ~ eV) - f(u>)) (11) 
J^ = 4vr 2 4 F |E J du\G^ F M\ 2 PNN^)p 2 LL-^)[f^-eV)-f{u J + eV)} (12) 



Ji corresponds to normal electron tunneling between electrodes, I2 is a net transfer of 
single electron with creation or annihilation of pairs as an intermediate state. -Z3 corresponds 
to a process in which electron from normal electron is converted to a hole in superconductor 
- branch crossing process in language of BTK theroy 2 , while I a is the Andreev tunneling. 

The differential conductance G(eV) = dI/d(eV) as a function of eV = Unm — Use is 
shown in the Fig. |HJ Clearly, if there is a spontaneous current in the ground state, the 
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FIG. 8: The total differential conductance for the solution with and without the spontaneous 
current. 



conductance peak is split, similarly as in the DOS. We could expect such behavior because 
G(eV) is proportional to the DOS at the Fermi energy. And again this effect could be 
observable in the tunneling experiments. 

We have also extracted Andreev conductance form the total one and ploted in the Fig. 
El We can see that conductance associated with the Andreev processes is strongly enhanced 
when the current flows in the ground state. Unfortunately it could be very difficult exper- 
imentally measure Andreev conductance only. Despite the fact that for energies less than 
SC gap the only allowed process is Andreev reflection, as in the point contact geometry, it 
doesn't work in our system. The problem is that even at very low energies there is a finite 
DOS at the Fermi level due to ferromagnet. Naturally the pairing amplitude is induced in 
FM slab but this is not true energy gap in the quasiparticle spectra and we always deal 
with some single electron processes in tunneling events. 
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FIG. 9: Corresponding Andreev differential conductance for the solution with and without the 
spontaneous current. 

VI. 2D FFLO STATE 



Before closing discussion on the spontaneous current we wish to make a remark regarding 
the nature of the ground state in our system. To begin with we recall that recently it has 
been predicted 36] that under certain conditions a 3D-FFLO state is energetically more 
favorable than usual ID state. The 3D state manifests itself in oscillatory behavior of 
the pairing amplitude not only in the direction perpendicular to the interface but also in 
direction parallel to it. Moreover, changing the thickness of the FM slab, one can switch 
the ground state of the system between 3D and 1D-FFLO state 0, 

The current carrying ground state of our system can be interpreted as a 2D-FFLO 
state. The argument is as follows: The oscillations of the pairing amplitude in the direction 
perpendicular to the interface occur regardless whether the spontaneous current flows or not. 
Within the FFLO theory j?l ls| , the period of the oscillations is related to the x-component 
of the center of mass momentum of the Cooper pair Q. On the FM side of our model the 
FFLO periodicity is governed by Q = (2E ex /vF)^f L , where vp is the Fermi velocity vector. 
This can be interpreted as the usual 1D-FFLO state in confined geometry. On the other 
hand, when the current flows parallel to the interface, there is a finite vector potential in the 
y-direction. This can be regarded as a y-component of the Q-vector. So one can say that 
when the spontaneous current flows, the 2D-FFLO state is realized. Moreover when the 
FM thickness is changed the ground state of the system is switched between 2D- and ID- 



state, which manifests itself in spontaneous current flow or in the lack of it. In the present 
calculations this vector was found during the self-consistency procedure, as it is related to 
the vector potential in the y-direction. Moreover, the effective Q y changes its value from 
layer to layer leading to inhomogeneous FFLO-like state in both dimensions. 

VII. CONCLUSIONS 

The competition between ferromagnetism and superconductivity in FM/SC heterostruc- 
tures give raise to the Fulde - Ferrell - Larkin - Ovchinnikov (FFLO) state in these systems. 
The original bulk FFLO state manifests itself in a spatial oscillations of the SC order 
parameter as well as in spontaneously generated currents flowing in the ground state of 
the system. We have argued that a very interesting version of this phenomenon accures in 
FM/SC proximity systems. In short, due to the proximity effect and Andreev reflections 
at the FM/SC interface, the Andreev bound states appear in the quasiparticle spectrum. 
These states can be shifted to the zero energy by tuning the exchange splitting or the thick- 
ness of the ferro magnet, thus they became zero-energy mid-gap states which lead to various 
interesting effects. In particular, the occurence of spontaneous currents in the ground state 
can be related to the zero-energy states, as in the case of high-T c superconductors. It seems 
that some combination of both phenomena is realized in a real systems. The fact that 
oscillatory behavior of SC order parameter is strongly correlated with the crossing of the 
Andreev bound states through Fermi energy and the generation of the spontaneous currents 
further support FFLO - Andreev bound states picture. The experimental confirmation of 
the existence of the spontaneous (spin polarized) currents in the ground state would support 
the FFLO - Andreev bound states scenario in these structures. 
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